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ABSTRACT 

We argue that the ‘changing look’ AGN recently reported by LaMassa et al. could be a lumi¬ 
nous flare produced by the tidal disruption of a super-solar mass star passing just a few gravita¬ 
tional radii outside the event horizon of a ^ 10 ®Mq nuclear black hole. This flare occurred in 
a massive, star forming galaxy at redshift z = 0.312, robustly characterized thanks to repeated 
late-time photometric and spectroscopic observations. By taking difference-photometry of the 
well sampled multi-year SDSS Stripe-82 light-curve, we are able to probe the evolution of the 
nuclear spectrum over the course of the outburst. The tidal disruption event (TDE) interpre¬ 
tation is consistent with the very rapid rise and the decay time of the flare, which displays 
an evolution consistent with the well-known behaviour (with a clear superimposed re¬ 
brightening flare). Our analysis places constraints on the physical properties of the TDE, such 
as the putative disrupted star’s mass and orbital parameters, as well as the size and tempera¬ 
ture of the emitting material. The properties of the broad and narrow emission lines observed 
in two epochs of SDSS spectra provide further constraints on the circum-nuclear structure, 
and could be indicative that the system hosted a moderate-luminosity AGN as recently as a 
few 10^ years ago, and is likely undergoing residual accretion as late as ten years after peak, 
as seen from the broad Ha emission line. We discuss the complex interplay between tidal 
disruption events and gas accretion episodes in galactic nuclei, highlighting the implications 
for future TDE searches and for estimates of their intrinsic rates. 
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1 INTRODUCTION 


It is widely accepted that super-massive black holes (SMBH) in 
the nuclei of galaxies grew over cosmolo gical times mainly by 
accre ting matter from their surroundings JsoltanI 1 19821 : iMerlonil 
I2OI5L and references therein). Decades of wide and deep surveys 
at different wavelengths have endowed us with a robust view of 
this growth process by constraining the luminosity function of Ac¬ 
tive Gal actic Nuclei (AGN) over most of the age of the Universe 
(see e.g. Hasingeret^ j|200^;|Ho;DGnsetalj20oi7yAirf e^ lj20ld: 


Ueda et alT 2014 : Buchner et alj 201^ Merloni T2OI5I : Aird et alj 

2OI5I : iBrandt & Alexander! 1201 Sl . and references therein). More¬ 


over, by studying simultaneously and coherently the properties of 
the AGN host galaxies within multi-wavelength surveys, we can 
now infer the overall duty-cycle of the AGN phenomenon, i.e. 
measure the fraction of time a typical galaxy spends in an ac¬ 
tive phase above a given (nuclear) luminosity jAird et alj|2012l : 




ongiomo 


et alJ20U t iHickox et alj20i^ . However, the results of 


these population-based (‘snapshot’) surveys do not provide any in- 
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formation on the behaviour of accreting black holes on timescales 
smaller than those over which galaxies evolve (~ billions of years). 
This is due to our uncertain knowledge of the fuelling mechanisms 
of AGN: for every galaxy of given properties, we still do not know 
whether black hole growth is rare and long-lived, or frequent and 
short-lived. 


In fact, understanding how black holes get their fuel at rates 
sufficient to power the observed AGN population, despite the huge 
angular momentum barriers present, is a long-standing goa l in the 
field of galaxy dynamics (see e.g. IShlosman et al.l 119901 : llogeel 
1 20061 : iHopkins & Ouata^l20ld , and reference therein). As sum¬ 
marized in the recent review by I Alexander & HickoxI 1 20121) . the 
current consensus is that gas inflow from kpc scales down to the 
central ~ 100 pc region occurs in all gas-rich galaxies, with the 
main driving mechanisms being a variety of externally-excited 
(mergers, fly-by) or internally-triggered (secular) gravitational in¬ 
stabilities, which help to shed a large fraction of the angular mo¬ 
mentum of the gas. On the other hand, on scales of order 1-10 pc 
(i.e. close to, or within, the gravitational sphere of influence of 
the central black hole), these large-scale instabilities are less and 
less effective. Therefore, other dynamical mechanisms must be in- 
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voked, such as dynamical friction, tidal disruption of clouds or star 
clusters, bars-within-bars, self-gravitating discs, eccentric discs or 
singled-armed spiral modes, driven by t he complex, time-varyin 
potential of the black hole, stars and gas i Jogee | 200d : lBekkill200l 


IShlosman et al ]|l990l:lHoDkins & OuataerdbOlOir 

Such complexity defies any simple modelization, as is appar¬ 
ent in the implementations of black hole growth within theoretical 
models of structure formation. There, lacking enough spatial reso¬ 
lution to model ab initio gas and stellar dynamics within galaxies, 
a number of different mechanisms have been postulated as primary 
drivers of nuclear accretion, often with equally satisfactory results 
(see e .g. lHopkins et alJl200^ : lHirschmann et alJl2012l : lMenci et all 
[Bn for recent treatments). 

Some observational attempts to infer recurrence times and du¬ 
rations of AGN accretion episodes have been made from indirect 
(typically space-resolved) tracers of past activity. A classic example 
is the light echo of Sgr A* that is apparently reveal ed by refection 
from the inner molecul ar zone of our own galaxy jSunvaev et alj 
ll993l : |Ponti et al.ll2013tl . Still within the Galaxy, excess Ha emis¬ 
sion along the Magellanic Stream has been claimed to be the result 
of illumination from Sgr A* about 1 — 3 x 10® years ago, when 
the black h ole was active at a level of 0.0 3-0.3 times its Eddington 
luminosity teland-Hawthorn et alj|201^ . Arguments for powerful 
past AGN events in nearby galaxies have been based on studies 
of extended emission-line regions, which can be used to trace the 
history of AGN emission over timescales of the order of the light 
travel ti me from the nucleus to the gas (typically 10'*-10® years , 
see, e.g., Dadina et 'Zl l20ld: iKeel et alJlBo: iGagne et al.ll20l4 
IPavies et al.l l2015h . Radio galaxies offer alternative routes into 
AGN variability, via the analysis of the morphology of the large- 
scale radio emission. Various pieces of indirect evidence of inter- 
mittency have been presented, such as the ripples and shock waves 
detected in the X-r ay emitting atmosphere around Virgo A/M87 
dporman et alj|20 0^. or the number versus size counts of sma ll ra¬ 


dio galaxies ( Reynolds & Begelmanlll997l : ICzemv et alj|200d) . Fi¬ 


nally, the decrease in the optical depth for Lya photons due to inter¬ 
vening absorbers along line s of sight to high redshift quasars (the 
so-called ‘proximity effect’ [Carswell etBll 19871 : ILu & Yull201lh 
can also be used to set constraints on the duration of individual 
QSO episodes. 

Yet another complicating factor is that at least some part of the 
fuel supply for AGN must come from stars, which abound within 


orbits passing close enough to the central black hole 

Hills 

1975; 

Gurzadian & Ozemoil 198ll: ICarter & Luminedll982l: 

Rees 

1988; 

Milosavlievic et alj 20061). The frequency of such eve 

nts depends 
nd is a non- 

ank & Rees 

on the stellar dynamical properties of galactic nuclei, a 
trivial outcome of a series of complex processes JFi 

1976l;lMaeorrian & Tremainell999l:lAlexandeil2012l: MerrittI20131: 


Vasilievi 120141) . Indeed, there is substantial uncertainty in the 
true rate of tidal stellar disruption in galactic nuclei. On the 
one hand, recent theoretical models, which account for realis¬ 
tic stellar dynamical models in galactic nuclei of different sizes 
and masses, have converged towards rates as high as Etdb ~ 
10“'^ yr~^, with j ust a weak dependence on the black hole (o r 
host galaxy) mass dWang & Merrittll2004l : IStone & Metzgeill2014) . 
On the other hand, observational studies have mostly reported 
constraints on the TDE rate that are up to one order of mag¬ 
nitud e lower (see, e.g., Donlwet ah! I2OO2I : Ivan Velzen & FarrM 
I2OI4I : iKhabibullin & Sazonovll2014l) . This apparent contradiction 
should not be surprising, as, observationally, the field is still in 
its infancy; only around two dozen TDE candidates have been 


identified so far by mean s of X-ray dKomoss a & Bade l 1 19991: 
Komossa & Greinei 1999al: Esquei et alJ I2OO8I: ICappelluti et alj 


200^ MaksymeUdJgOl^ Bloom et^l j2()H : Burrows^etalj2()lT|: 


Saxton et al 2()lj| Jcenkort ^ 2012lJ : Khabibullin & Sazonov 

201A ) ■ UV iGezari et al.j 200^ 200S) and optical 1 van^elzenet^ 


201T jGezarirt alfeoddA rcavi et al ]|2014l:IChornock et alfcOll 


Holoien et alj 2014 ) observations (see Gezaril 2014 , for a recent 
overview). In particular, the full extent of the various selection 
effects that may plague TDE selection in different bands of the 
electromagnetic spectrum have yet to be properly characterized. 
Nevertheless, the rate of discovery is increasing dramatically with 
the advent of large wide-area optical time-domain surveys (mostly 
driven by supernovae sea rches), such as the Catalina Real-Time 
Transient Survey (CRTS, DrAe et ^ l2009l) . the Palomar Tran¬ 
sient Factory (PTE, iRau et alj|2009l) . PanSTARRS dKaiser et alj 
bold) and the All-Sky A utomated Survey for Supemovae (ASA- 
SN. Ishappee et aLll2014l) . The future promises even more rapid 
advances with next generation wide area X -ray (SRG/eROSITA, 
iMerloni e t alj|2012l: ^abibullin et alj|2014l) and optical ( Skymap- 


per. 


Keller et alj2007l. ZTitH and LSST.Ilvezic et alj2008l) 


)surveys. 


In this work, we focus on just one particular example: a 
2 = 0.312 galaxy within the Sloan Digital Sky Survey (SDSS) 
‘Stripe-82’ are^^ which has undergone a dramatic and rapid change 
in its nuclear emission, from type-1 (broad line) AGN-like, to type- 
1.9 (sho wing only a weak broad Ha emission line) within a 10 year 
interval. iLaMassa et alj d2015l) . who first reported this serendipi- 
tous discovery, rightly point out that suc h a dramatic change is rare 
jPadina et alj201(il : lDennev et alj20l^ . and can provide important 
clues on the nature of black hole fuelling in galactic nuclei. Here we 
perform an in-depth analysis of the available (mostly public) data 
on this source (section and O, and we argue, based on the out¬ 
come of our analysis, that its overall light-curve in the last decades 
is consistent with the tidal disruption of a main-sequence star (with 
mass M* > IMq) by a « 10® Mq black hole hosted in the nu¬ 
cleus of a massive, star-forming galaxy (section O. Our conclu¬ 
sions are mainly based on the properties of the optical light-curve, 
and the overall energetics of the event. However, the full picture 
reveals the complexity of this galactic nucleus, as we demonstrate 
that (i) the material giving rise to the broad emission lines illumi¬ 
nated by the observed flare cannot be a distant part of the stellar 
debris (section ISTTT ). and (ii) the narrow emission lines could indi¬ 
cate that about 10"^ years ago the central black hole was active at 
a level similar to that observed (section lY2l ). This prompted us, in 
section l53] to develop a simple unified model for AGN fuelling, in 
order to estimate the relative occurrence of gaseous accretion and 
tidal disruption flares in nearby galactic nuclei. Finally, we draw 
our conclusions in section 


2 DATA ANALYSIS 

2.1 Source identification and long term evolution 

The main focus of this paper is the object 
SDSS J015957.64-H003310.5, hereafter SDSS J0159-H0033. 
SDSS J0159-I-0033 first came to our attention (before being 
independently reported by ILaMassa et aklboisl) whilst creating a 


^ http://www.ptf.caltech.edu/ztf 

^ Stripe-82 is a 109 X 2.5 degree survey field that lies along the celestial 
Equator within the Southern Galactic Cap, which has been imaged multiple 
times by the SDSS camera 
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reference sample of bright X-ray sources to guide target selection 
for the SDSS-IV/SPIDERS survey. (Merloni et al., in prep., 
see also www.sdss.org/surveys/). Starting from a parent 
sample of bright {Fo, 2 - 2 keV > 10“^® erg cm“^ s“^) X-ray 

sources taken from the 3XMM serendipitous source catalog we 
found 557 objects having optica l spectroscopic ob servations in the 
latest SDSS data release (DR12. lAlam et alj2015h . A small subset 
(35/557) of these had a good quality spectrum available from 
both the SDSS- I/II phase (using the original SDSS spectrograph, 
lYork et alfcOOOh . as well as a spectrum obtained durin g the SPSS- 
Ill ph ase (using the upgraded BOSS spectrograph, Smee et al] 
l2013h . Within this subset of 35 objects there was a single one, 
SDSS J0159-I-0033, for which the automated classification derived 
from the SDSS-I/II spectrum (CLASS='QSO' ) disagreed with 
that derived from the SDSS-III spectrum (CLASS='GALAXY'). 
SDSS J0I59-I-0033 appeared in our original X-ray sample due 
to its serendipitous detection in an X MM-Newton observati on of 
the nearby luminous QSO Mkn 1014 dMarkarian et al.|[T9^ . The 
apparent disparity between the two spectroscopic classifications of 
this object prompted further investigation. 

SDSS J0159-I-0033 was first imaged by the SDSS cam¬ 
era on 1998 September 25, and was subsequently selected 
for spectroscopic follow up within the low redshift ugri 
colour selected QSO targ et category (specifically with target 
bit T ARGET_QSO-SKIRT, iRichards et'Zll2002l : IStoughton et alJ 
l2002h . The pipeline measurements derived from the SDSS photom¬ 
etry found SDSS J0159-I-0033 to be a marginally resolved galaxy 
having de-reddened ugri colours (u — g = 0.53, g — r = 0.79, 
r — i = 0.31) that placed it outside the main stellar locus, blue- 
ward of the colour regions used to reject ‘normal’ galaxies, and 
also outside additional colour boxes designed to reject rarer stellar 
contaminants. The optical spectrum taken on 2000 November 23 
(MID 51871) is classified as a 2 = 0.312 broad line QSO, due the 
unmistakable presence of a blue continuum and broad Balmer lines 
(Ha and H/3, see section lZ^ for more details). 

On 2010 January 5 (MID 55201), SDSS J0159-t0033 was ob¬ 
served spect roscopically by the B aryon Acoustic Oscillations Sur¬ 
vey (BOSS. iDawson et alJlioTjh within the SDSS-III project, as 
part of a small random sub-set of known broad line AGN from 
SDSS that were re-observed in order to define QSO spectral tem¬ 
plates for BOSS (targeting flag TEMPLATE_QSO_SDSS). However, 
the resulting spectrum revealed a 2 = 0.312 star-forming galaxy 
with essentially no blue continuum, and no prominent broad emis¬ 
sion lines, apart from a weak broad shoulder to the Ha line. The 
BOSS pipeline measured a velocity dispersion of cr = 169 ± 24 
km/s. 

Five X-ray observations of the field containing 
SDSS JO 159-1-0033 exist in public archives. The oldest of 
these is a Roentgen Satellite (ROSAT) All-Sky-Survey scan 
from 1991 January, which failed to detect the source, placing 
a 3cr upper limit of 3.2 x 10“^cnts“^, corresponding to an 
unabsorbed (i.e. c orrected for the Gal actic column density of 
2.3 X 10^° cm~^. iKalberla et alj [20051) X-ray flux limit in the 
0.1—2.4 keV band Tb.i— 2.4kev < 8 x 10““ergcm“^s"^ (fora 
Tx = 2.1 power-law spectrum). Two deeper (6.1 and 1.6 ksec), 
ROSAT-PSPC pointed observations of nearby sources (Mkn 1014 
and MS0158.5-H0019) on 1992 January 17 and 1992 July 24 also 
covered the location of SDSS J0159-I-0033, placing more stringent 
upper limits on its flux of To.i-2.4keV < 2.6 X 10 and 


^ http://xminssc-www.star.le.ac.uk 


Tc).i- 2 . 4 keV < 4.6 X 10“^® ergcm“^ s“^ respectively. In contrast, 
the XMM-Newton observations of the same region performed 
in 2000 July (aimed at the nearby luminous QSO Mrk 1014), 
caught SDSS J0159-t0033 in an X-ray bright state, with a flux 
of Ti-iokeV = 2.6 X 10“^^ ergcm“^ s“^, and a spectral shape 
typical of accreting super-massive black holes, characterized by 
a power-law with slope Tx = 2.1, and no sign of significan t 
absorption at the redshift of the sourc^ jLaMassa et alj |20T^ . 
Finally, a Chandra observation of Mkn 1014 serendipitously 
covered SDSS J0159-I-0033 in 2005, and measured it to be a factor 
of ~7 fainter than in 2000, but again with a power-law spectrum 
consistent with no absorption at the redshift of the source. 

To place the observed high-amplitude variability, detected in 
both X-rays and SDSS optical spectra, into a long-term context, we 
have collected publicly available UV and optical photometric data- 
points of the source spanning the last ~30 years, which we display 
in Fig.[T] 

There is no evidence for significant photometric vari¬ 
ability after 2005. The spectral energy distribution (SED) of 
SDSS JO 159-1-0033 derived from the most recent Stripe-82 pho¬ 
tometric observations (autumn 2007) shows a great deal of con¬ 
sistency with the BOSS spectrum taken in 2010 (see Eig. 
and section 12.21 below). As a further check, we obtained 7- 
band (g’r’i’z’JHK) GROND jGreiner et al.ll200^ observations of 
SDSS J0159-H0033 on 2014 November 04, and found that the SED 
was still dominated by the host galaxy, fully consistent with all pho¬ 
tometric observations since about 2005, including public photome¬ 
try from wide-field, high-cadence, optical transient surveys (CTRS 
and PTE), and als o confirmed by the P alomar spectrum taken in 
2014 (presented bv iLaMassa et ai.lEoi.5h . 

More interestingly, we find that also all available data-points 
before 1998 are consistent with the optical emission of the source 
being dominated by its host galaxy, albeit with larger uncertain¬ 
ties due to the poorer sensitivity and calibration of photographic 
plates. The older photometric measurements shown in fig. [T] were 
derived from th e Palomar and UK- Schmidt photographic plate sur¬ 
veys (GSC2.3.2. lLasker et al.l2008h . and have been converted to the 
nearest SDSS filters according to the procedure described in Ap¬ 
pendix \E These also include a 1983 V-band observation (’’Palomar 
Quick-V Northern Survey”, not shown in Fig. [T} which translates 
to r’= 19.17 ± 0.36 (see Tab. lAll in Appendix IaI. Uncertainties 
on those fluxes, however, are such that it is not possible to exclude 
that SDSS J0159-I-0033 was in a state of (low-level) AGN activity 
even before the flare began. 

In summary, the only period during the last ~30 years when 
SDSS 10159-1-0033 was significantly brighter than its late-time 
(post-2005) constant level (both in X-rays and optical bands) was 
during the period between 1998 and 2005, suggesting a one-off, 
high amplitude, flaring episode, rather than a more typical AGN 
variabi lity pattern, often mod elled as red-noise or ‘‘damped random 
walk” dMacLeod et alJl^loL and references therein). 

2.2 Spectral analysis 

lUaMassa et alj ( 1201 5h have presented a detailed analysis of the 
2000 and 2010 spectra, and derived the mass of the central 
black hole under the assumption that the observed broad Ha and 
H/I emission lines are produced by virialized material in orbit 

^ Within this model, this corresponds to Uq i_ 2 . 4 keV = 6.4 X 10“^^ 
ergcm“^ s“^. 
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MJD (days) 

Figure 1. Top Panel: Long term optical light-curve for SDSS J0159+0033. The SDSS Stripe-82 photometric measurements in the u’g’r’i’z’ filters are shown 
with small open circle symbols (blue, green, red, magenta and black respectively), with the yearly medians for each filter connected by solid lines. Stripe- 
82 epochs not meeting the data-quality criteria described in the text are omitted. The spectro-photometric measurements derived from the two epochs of 
SDSS spectroscopy (also for the iPg’r’i’z’ filters) are shown with large open boxes (offset by —0.3 magnitudes to align roughly with the contemporary 
SDSS photometric measurements). Ol der photometric measurements derived from the Palomar and UK-Schmidt photographic plate surveys, are individually 
labelled fGSC2.3.2 lLasker et al.l2008h . They have been converted to the nearest SDSS filters according to the procedure described in Appendix m The yearly 
median and RMS of the unfiltered CRTS (Catalina T ransient Survey, Drake et al .l2009h measurements are shown with (grey) filled pentagons. The median and 
RMS of the public RTF (Palomar T ransient Factory, | Rau et al.ll20Q9 ) data in the g' and R filters are shown with downward-pointing triangles (green and red 
respectiv ely). GROND phot ometry ^Greiner et al.l2008l) in the e’r’i'z' bands is shown with filled square symbols (green, red, magenta and black respectively). 
GALEX ^Martin et al.ll2005h data-points are shown with filled and open light blue diamonds for the NUV and FUV bands respectively. The epochs of the 
ROSAT, XMM-Newton and Chandra X-ray observations are shown with vertical dashed lines. The Bottom panel: shows the X-ray lightcurve. 


around the central black hole. Adopting the relation between ob- 
served line width, cont i nuum lum inosity and b l ack h ole mass of 
IVestergaard & PetersonI < l2006h and iGreene et alj feOld) for H/3 and 
Ha, respectively, they derived a mass of about 1.7 X 10®M©. Here 
we briefly describe the results of our independent, detailed, spectral 
analysis. 

The SDSS and BOSS spectra were obtained through different 
sized apertures (3.0 and 2.0 arcsec diameter respectively). There¬ 
fore we first corrected the spectro-photometric calibration of the 
SDSS (2000) and BOSS (2010) spectra by matching the i’-band 
spectro-photometric flux measurements (which, at the redshift of 
the source, z =0.312, is devoid of strong emission lines), to the 
nearest (in time) Stripe-82 i’-band photometric measurements. For 
the SDSS 2000 spectrum this was the single Stripe-82 data-point 
closest to the peak of the outburst, and for the BOSS 2010 spec¬ 
trum we took a weighted average over the autumn 2007 i’-band 
photometric measurements (see figure[2]l. 

We then fitted the BOSS 2010 spectrum in the regions of the 
most prominent emission lines (Oil, 3650A < A < 3800A ; H/3, 
4700A < a < SIOOA; Ha, 6480A < A < 6750A). For the host 


galaxy continuum and absorption lines, we used the SSP (single 
stell ar population) high-resolution e volutionary model templates 
from lOonzalez Delgado et ^ ( l2005h . the preferred model (from a 
comparison with the broad band SED) being a Z = 0.019 (i.e. 
~ solar metallicity) ‘Padova’ isochrone template with an age of 
2.5 Gyr. The normalization of this component corresponds to a to¬ 
tal mass of formed stars of « 1.1 x lO^^M©, which, taking into 
account stellar mass loss, would imply a stellar mass of the host 
galaxy of about « 8 x lO^^MeQ. The corresponding model SED 
is shown in Eig.|2 and in the bottom panels of Eig. |4] To the stel¬ 
lar component, we added a power-law continuum to represent any 

® We note that we are not aiming at a precise determination of the age of the 
stellar population. Indeed, if we chose a younger SSP from the same library 
(for example a 1.6 Gyr old one, in agreement with LaMassa et al. 2015), we 
would recover an equally good fit to the nan'ow wavelength ranges around 
the emission lines we are interested in, with all emission line parameters 
within the errors, and a stellar mass of the host of 7 x 10^®M©. This 
shows that, on the other hand, the total stellar mass of the system is quite 
robustly determined by the SSP normalization. 
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Figure 2. The two dark grey lines show the calibrated spectra (\Lx, in 
units of 10"^^ erg/s) of SDSS J0159+0033 from 2000 (upper curve) and 
from 2010 (lower curve). Superimposed ai'e the Stripe-82 u’g'r’i'z' photo¬ 
metric data points from near the outburst peak (2000, black solid circles), 
and the weighted average of the last season of Stripe-82 photometry (2007, 
black solid stars). The red solid line is the best-lit model to the continuum 
emission of the 2010 spectrum, which as described in detail in the text, is 
the sum of a SSP host galaxy spectrum (solid puiple line) plus a power- 
law continuum, whose slope has been fixed at Fiambda oc (purple 

dashed line). Finally, the open star symbols show the implied ‘pure’ galaxy 
SED that we use to study the nuclear fiare evolution. 


residual nuclear emission, which is required to fit the blue end of 
the observed spectrum. However, the slope of such a component 
is barely constrained by the BOSS 2010 spectrum, and so we de¬ 
cided to_fixJt_to_that_e?™cte^or the ‘canonical’ viscous accretion 
disc JShakura & Sunva^l 19731) a = —1/3, where Fx oc 
Changing the slope of the power-law component within the range 
observed from the flare emission during the outburst (see section[3] 
below) changes our results only marginally, and not qualitatively. 

We fit the emission lines with the MPFIT IDL routine 
jMarkwardtl 1 20091) . using a set of Gaussians, as described below; 
the best-fit parameters for the emission lines (of both the SDSS 
2000 and BOSS 2010 spectra) are shown in Table[T] 

For the Fl/3 region in the BOSS 2010 spectrum (see top 
left plot in fig. S, one narrow component is sufficient to fit the 
Balmer line, while the [Olll] complex (where we fixed the flux 
of [O III] 4959 A to be 1 /3 of the [O III] 5007 A flux) requires 
an additional, blue-shifted broad component, typically associated 
with outflows in the narrow-line region ionized gas. This blue wing 
has a measured FWFIM of 774±87 kms“^ and is blue-shifted 
by about 100±50 kms“^ with respect to the narrow component. 
Such a width for the broad component of [Olll] is well within 
the average for Seyfer t 2 galaxies in SDSS of similar luminosity 
as SDSS J0159-I-0033 jMullanev et al.ll^l3h . but not so high as to 


place it squarely among the secure sources of AGN-dri ven outflows 
(typic ally requiring FWHM > 1000 kms“^, see e.g. iBrusa et alj 
[ 2 ^ . 

In the Ha region of the BOSS 2010 spectrum (see top right 
plot in fig. 0, in addition to the narrow Ha line and [N II] complex 
(where, again, we fixed the flux of [N II] 6548 A to be 1/3 of the 
[N II] 6583 A flux ), a ver y broad line is required, as already noticed 
bv iLaMassa et al.l l l2015tl . 

The narrow lines observed in the BOSS 2010 spectrum can 
be used to investigate the source of ionizing radiation that excited 
them. We plot in Fig. [3] the BPT ( [Baldwin et alJligsU) diagnostic 
diagram for our measurements of the narrow emission lines ([O III] 
H/3, [NII] and Ha) as well as those derived from the analysis of the 
three spectra presented bv iLaMassa et alj ( 12015l) (the 2000 SDSS 
and 2010 BOSS spectra, as well as a later Palomar DBSP spec¬ 
trum taken in 2014). As expected, since we fit two components to 
th e [Olll] lines, we find different H/3/[0 III] flux ratios to those 
of iLaMassa et al.l ( l2015h . who adopt a single component model, 
whilst the [N Il]/Ha and [S Il]/Ha line ratio measurements are in 
agreement. If we only consider the narrow components of [Olll], 
SDSS J0159-(-0033 appears to lie well within the ‘transition’ re¬ 
gion between star-forming galaxies and AGN. On the other hand, 
if we sum together both narrow and broad components of the [O III] 
line, the object moves into the AGN-dominated part of the diagram. 
Thus, it appears that the line emitting material was probably ionized 
by a source harder than that associated with pure star-formation. 
AGN activity would of course provide a sufficiently hard ionizing 
source, but the exact amount of any putative AGN contribution to 
the observed narrow line emission is harder to disentangle. We dis¬ 
cuss the consequence of the narrow emission lines properties for 
the interpretation of the flare in section lS^ 

Moreover, the strong [Oil] emission line can also be used to 
infer an approximat e star formation rate (SF R) of the host. We fol¬ 
low the procedure of isilverman et alj ( l2009h to correct the [Oil] for 
any possible AGN component, by assuming that all observed [OIII] 
emission is AGN-driven, thus obtaining a lower limit to the derived 
SFR. We assume a fixed [OII]/[OIII]= 0.21 ratio for AGN-only 
excited emission lines, and use the difference I/fom — I/[oii],agn 
to derive a SFR (see eq. (2) of I Silverman et alj|2009t) . Under this 
assumption, we obtain log SFR« 1.6(MQ/yr). 

A direct comparison between the emission lines seen in the 
2010 BOSS and 2000 SDSS spectra can reveal the variable compo¬ 
nents, and help determine their physical origin. In our analysis of 
the 2000 SDSS spectrum, taken very close in time to the photomet¬ 
ric peak of the light-curve (see Fig.[TJ, we have first assumed that 
the host galaxy continuum and the narrow emission lines seen in 
the 2010 BOSS spectrum have not changed significantly over the 
intervening ten years. This baseline constant spectrum is shown as 
a dashed red line in the middle panels of fig. (4] We then take the 
difference of the two spectra (difference = 2000 SDSS spectrum 
- 2010 host galaxy emission model), plotted in the bottom pan¬ 
els of Fig. [4] In the H/3 region, a clear rising continuum is visible, 
with a prominent broad H/3 emission line. We fix the power-law 
continuum slope to that measured over the entire spectral range in 
the difference spectrum, a = —0.5, leaving its normalization free, 
and measure the flux and FWHM of the broad lines (also reported 
in table [TJ. The [Olll] emission region is noisy, but the spectrum 
taken at the peak of the outburst does also show a broad compo¬ 
nent, which is consistent with having the same properties of that 
observed ten years later (FWHM of 753±80 km/s). 

We can use the measured parameters of the broad emission 
lines to estimate the mass of the central black hole, using the stan- 
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Figure 3. Line emission diagnostic fBPT. lBaldwin et al.ll98lh diagrams for 
the three spectroscopic observations (2000 SDSS, black circ le; 2010 BOSS, 
triangl e and 2014, square) of SDSS J0159+0033 presented in lLaMassa et alj 
i2015l) . The red triangles show the results of our own analysis of the 2010 
spectrum: the filled symbol corresponds to the case in which we only con¬ 
sider the nan'ow [O III] line component, the empty one for the sum of broad 
plus narro w [O III] line compon ents. In each panel the dot-dashed line rep¬ 
resents the lKewlev et alj <20061) separation between pure AGN and the rest 
of the galaxy population. In the left panel the dashed line is the separation 
between star-forming galaxies and ‘transition’ objects (i.e. those with ion- 
izing photons coming fr om both star-forming regions and AGN) defined by 
iKauffmann et aDES). 


dard ‘single epoch’ virial method (se e e.g. Peterson et al. [200^ 
and references therein). We adopt the I Greene & Hoi \ 2007 1 scal¬ 
ing between black hole mass, Ha FWHM and Ha line luminos¬ 
ity to derive log(MBH/A/ 0 ) ~ 7.8 and /Mq) ~ 8.2 

for the 2000 and 2010 spectra, respectively. Scaling laws to ob¬ 
tain black hole masses from the measured width of the H/3 line 
and the continuum luminosity at 5100 A (in 2000 measured to be 
AI/5100 = 6.8 X 10^^ erg s~^) h ave been published by many au - 
thors. Adopting the calibration of IVestergaard & Petersoii| j[200^ 
we de rive log(MBH/MQ) ~ 8.1, while following Greene et alj 
LOlOh we get log(MBH/Af 0 ) ~ 8.2. All these estimates come 
with both statistical and systematic uncertainties of about 0.3 dex, 
due to the uncertain calibration and the unknown ge ometry of the 
broad line region itself (see e.g. IShen & Kellvll20la) . Encouraged 
by the c onsistency among all the virial black hole mass estimators 
(see also lLaMassa et alj2015tl . we consider these estimates robust, 
and, for the remaining of the paper, assume a fiducial SMBH mass 


of 10 ® Mq. 

In sections ISTI and we discuss the observed properties of 
the emission line regions within the context of our interpretation of 
the observed flare from SDSS J0159+0033. 


3 THE BOLOMETRIC LIGHT-CURVE OF THE FLARE 

A change in the optical spectral properties of AGN as dramatic as 
that observed in SDSS J0159+0033, and taking place within just 
ten years, is very unusual. iLaMassa et alj ( 120151) have discussed a 
possible interpretation of these rapid changes of flux and spectral 
properties within the context of the extreme end of normal AGN 
variability. Given the implied size of the broad line region, and the 
observed X-ray spectra (both XMM-Newton and Chandra can be 
fitted with un-absorbed power-laws), they conclude that an obscu¬ 
ration event, in the form of a large-scale cloud passing in front of 


the central source and of the broad line region, is not consistent 
with the timescale of the observed fast flare, nor with the spectral 
slope of the 2005 Chandra spectrum. It appears, then, that the in¬ 
trinsic emission of the AGN is fading rapidly, possibly due to a 
decrease of the accretion rate onto the central black hole. However, 
in the context of standard accretion-disc theory, it is not straightfor¬ 
ward to explain both the rapid fading in the period 2000-2005 and, 
in particular, the very fast brightening of the source between 1998 
and 2000 ; viscous timescales of optically thick and geometrically 
thin accretion discs as close as 10 gravitational radii from a 10 ® M 0 
black hole, with the luminosity observed close to its peak, are of 
the order of a few years (see section|5]below). Larger discs, needed 
if the observed optical light is produced from viscous dissipation 
within the disc itself, would evolve on even longer timescales. 

It appears, then, that a physical process is required that can 
give rise to rapid enhancement of the accretion rate, a process that 
delivers large quantities of matter very close to the central black 
hole, where viscous (or thermal) times are the shortest. Moreover, 
such a process must be be highly intermittent, because, as we dis¬ 
cussed before, the long term light-curve of SDSS J0159+0033 is 
essentially dominated by this one single flare. In the following, we 
describe why we believe that a Tidal Disruption Event (TDE) is a 
much more natural interpretation for this peculiar transient event. 

We first analyze in more detail the evolution of the flare, as 
traced by the Stripe-82 lightcurve. Our starting point is the Stripe- 
82 photometry for SDSS J0159+0033 and two neighbouring ref¬ 
erence stars (SDSS J020002.27+003250.6, a r « 18 magnitude 
star about 1.2 arcmin away, and SDSS J015955.54+003419.6, a 
r « 16.1 star about 1.26 arcmin away), retrieved from the SDSS 
DR7 data serveB In the few cases where multiple data-points were 
found with the same MJD, we took daily averages, resulting in 70 
independent epochs of photometry. We then cleaned the light curve 
by excluding all epochs with poor photometry, as deduced by com¬ 
paring the flux of each reference star in each epoch and each band 
with their long-term average, and discarding all epochs where the 
average flux of the two stars in any band is more than 2a off its 
mean. This leaves us with 50 high quality epochs of photometry. 

We then calculate the ‘baseline’ non-varying host galaxy lu¬ 
minosity in each of the five SDSS filters by assuming that the 
weighted average of the light collected in the last 10 epochs of 
Stripe-82 data (autumn 2007) is due to the large-scale galaxy stel¬ 
lar emission, plus a small contribution from a nuclear power-law 
(which we constrained by analysing the 2010 BOSS spectrum, see 
section [Tl] above). The SED of the baseline ‘pure’ (i.e. with the 
nuclear power-law component removed) galaxy emission is shown 
in Eig.l^as empty stars. 

We then derive the ‘nuclear’ luminosit 3 lj| in each SDSS filter 
by subtracting the baseline host flux from the observed flux at each 
Stripe-82 epoch (see Table column (2) to ( 6 ), where we report 
the source nuclear luminosities based on SDSS modelMag). Ei- 
nally, we compute the total optical luminosity of the nuclear flare 
Lopt.nuc = Lxd\ ~ Lx^iAXi, where the integral ex¬ 
tends over the entire SDSS bandpass, and the summation is over 
the u’g’r’i'z’ SDSS f ilters, each characterized by a bandwidth AAi 
dEukurita et inil996h . 

The top panel of Fig.[^shows the time evolution of the optical 


° http://cas.sdss.org/stripe82/en/tools/crossid/crossid.asp 
^ All luminosities have been computed from the measured fluxes adopting 
a ACDM cosmology with Sim = 0.286, S2 a = 0.714 and Hq = 69.6 
km/s/Mpc. 
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Table 1. Emission line parameters derived from from our spectral analysis (see section lT^ for details). 



A(A) 

SDSS 2000 
FWHM(kms-l) 

-6iine(10'‘°ergs l) 

A(A) 

BOSS 2010 
FWHM(kms-l) 

Iiiine(10'‘°ergs 1) 

[Oil] 

3726.79“ 

468“ 

38.7“ 

3726.79±0.07 

468±13 

38.7±3.0 

H/3narrow 

4860.8“ 

309“ 

18.0“ 

4860.8±0.1 

309±13 

18.0±1.8 

R^broad 

4865.8±1.4 

4493±234 

103.9±5.1 

- 

- 

- 

[OlII] narrow 

5005.7“ 

308“ 

30.2“ 

5005.5±0.1 

308±13 

30.2±2.2 

[OlII] broad 

5006.0±0.5 

753±80 

18.6±2.6 

5003.8±0.6 

774±87 

21.6±5.1 

Hanarrow 

6562.26“ 

262.2“ 

67.9“ 

6562.26±0.05 

262.2±5.1 

67.9±4.1 

R^broad 

6559.9±0.9 

3408±110 

329±11 

6582.4±2.0 

6167±280 

143.0±6.7 

[Nil] 

6582.7“ 

276“ 

29.5“ 

6582.7±0.1 

276±13 

29.5±2.9 

[SII] 

6715.8“ 

358“ 

24.1“ 

6715.8±0.2 

358±18 

24.1 ±2.2 


“Line pai'ameter fixed at the 2010 values. 
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Rest—Frame Wavelength [AA] 


H„, [Nil] 
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Figure 4. Left: Spectral analysis of the H/3 region. The top panel shows the calibrated BOSS 2010 spectrum with uncertainties (black histogram with grey error 
bars), fitted with a host galaxy continuum model from Gonzalez Delgado et at (2005) (a 2.5 Gyrs old Z = 0.019 SSP; top purple dashed line), a power-law 
continuum of slope q = —1/3 (bottom purple dashed line), nan'ow H/3 and [OlII] emission lines (green solid line) and a broad [OlII] component (blue 
solid line). The middle panel shows the observed, calibrated 2000 SDSS spectrum, together with the 2010 best fit host galaxy model continuum plus naiTow 
emission lines (red dashed line). The bottom panel shows the difference spectrum (2000-2010), where the rapidly time-varying emission components emerge. 
This difference spectrum is fitted with a power-law of slope a. = —0.5, where Fx oc (dashed purple line) plus broad H/3 and [O III] emission lines 

(cyan line). Right: Spectral analysis of the Ha region. The top panel shows the calibrated BOSS 2010 spectrum with uncertainties (black histogram with grey 
error bars), fitted with a host galaxy continuum model from Gonzalez Delgado et al. (2005) (a 2.5 Gyrs old Z = 0.019 SSP; puiple dashed line), a power-law 
continuum of slope a = —1/3 (bottom purple dashed line), a nan'ow Ha and [Nil] emission lines (green solid line) and a broad Ha component (blue solid 
line). The middle panel shows the observed, calibrated 2000 SDSS spectrum, together with the 2010 best fit model host galaxy continuum plus narrow emission 
lines (red dashed line). The bottom panel finally shows the difference spectrum (2000-2010), where the rapidly time-vaiying emission components emerge. 
This difference spectrum is fitted with a power-law continuum (dashed puiple line) plus broad Ha emission line (cyan line). 


nuclear luminosity I/opt,nuc, where each epoch is represented by 
a data point of a different colour. A simple power-law fit to the 
time evolution of the flare in its decay phase gives a slope of rioo = 
— 1.59dz0.05, only slightly shallower than the predictions of simple 
tidal disruption flare models (a power-law evolution in with slope 


—5/3 lReeslll9^: iLodato et al.ll2Q0^: iGuillochon & Ramirez-Rui3 

l2013h . 

Fig- m shows the evolution of the optical SED of the flare; for 
clarity, only the first 16 epochs of Stripe-82 photometry are shown, 
with the same colour code of Fig. [5] As a reference, the plot shows 
also the constant baseline host galaxy spectral energy distribution 
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Table 2. Nuclear (host-subtracted) light-curve of SDSS J0159-I-0033 in the u’g’r’i’z’ SDSS bands 


MJD 

erg/s) 

Lx,g>^Xgi 

(10'‘2 erg/s) 

^X,r' 

(lO^^erg/s) 

L\ 4/AA^/ 
(loi^ers/s) 

Lx,.'^X,, 

(lO^^erg/s) 

51081 

7.1± 1.0 

9.3±0.6 

4.0± 0.5 

1.3±0.5 

2.7± 1.5 

51819 

29.8± 1.3 

39.9± 0.8 

17.4± 0.6 

9.2± 0.5 

16.4± 1.6 

52225 

20.9± 1.1 

26.5±0.7 

10.5±0.5 

4.9± 0.5 

11.0± 1.4 

52288 

17.1± 1.1 

21.9± 0.7 

9.5± 0.6 

5.5±0.5 

11.5± 1.5 

52522 

13.7± 1.1 

14.0± 0.6 

7.1±0.5 

4.0± 0.5 

7.4± 1.6 

52551 

13.6± 1.1 

11.4± 0.6 

5.2± 0.5 

1.6± 0.5 

6.1± 1.6 

52558 

12.9± 1.0 

13.5±0.6 

5.3±0.5 

3.1±0.5 

-2.9± 2.5 

52576 

12.3± 1.0 

15.3±0.6 

5.1±0.5 

2.7± 0.5 

6.4± 1.6 

52577 

14.4± 1.0 

16.8±0.6 

6.9± 0.5 

4.4± 0.5 

9.5± 1.5 

52585 

12.2± 1.1 

15.6± 0.6 

5.4± 0.5 

2.9± 0.6 

7.3± 1.8 

52909 

18.4± 1.1 

21.6± 0.6 

9.0± 0.5 

4.3±0.5 

10.1± 1.5 

52910 

19.2± 1.1 

21.1±0.7 

9.8± 0.5 

4.4± 0.5 

9.6± 1.5 

52935 

15.9± 1.1 

17.6± 0.7 

7.4± 0.6 

3.1±0.6 

7.8± 1.7 

52963 

14.8± 1.0 

16.2± 0.6 

5.5± 0.5 

3.7± 0.5 

7.3± 1.4 

52971 

10.6± 1.1 

12.6± 0.6 

5.0± 0.5 

2.1±0.5 

4.0± 1.5 

53351 

9.1± 1.0 

8.0± 0.6 

3.9± 0.5 

2.1±0.6 

3.0± 1.7 

53634 

3.4± 3.0 

6.5± 1.7 

1.6± 1.1 

1.8± 0.7 

3.6± 1.8 

53655 

5.0± 1.1 

4.5± 0.7 

3.2± 0.6 

0.8± 0.5 

2.5± 1.7 

53665 

2.6± 2.1 

3.9± 1.0 

1.6± 0.7 

0.3±0.6 

1.5± 1.8 

53669 

5.1± 1.1 

3.9± 0.6 

1.7±0.5 

1.0± 0.5 

2.3± 1.6 

53671 

4.6± 1.4 

3.6± 0.7 

2.3± 0.7 

0.9± 0.6 

3.6± 1.7 

53676 

3.0± 1.2 

2.5± 0.6 

0.6± 0.6 

0.1±0.5 

1.3± 1.5 

53679 

2.3± 1.4 

1.9± 0.7 

-0.8± 0.7 

-0.3± 0.6 

-0.0± 1.8 

53700 

4.9± 1.3 

4.7± 0.7 

2.0± 0.6 

1.2± 0.6 

6.4± 1.9 

53705 

3.5± 1.3 

2.8± 0.7 

2.0± 0.6 

0.9± 0.6 

-0.8± 1.7 

53975 

3.1± 1.2 

4.5± 0.6 

2.2± 0.5 

1.5±0.5 

2.4± 1.7 

53994 

4.0± 1.2 

2.7± 0.6 

1.0± 0.5 

0.4± 0.5 

-3.0± 1.4 

54008 

2.8± 1.1 

3.0± 0.6 

1.8±0.5 

1.4± 0.5 

1.4± 1.7 

54010 

2.7± 0.9 

3.1±0.6 

-0.3± 1.0 

1.0± 0.5 

0.1± 1.4 

54012 

4.1± 1.8 

0.0± 1.1 

-1.6± 0.7 

-2.4± 0.6 

0.2± 1.7 

54020 

2.3± 1.7 

5.8± 1.0 

4.8± 0.7 

2.6± 0.7 

4.5±2.0 

54024 

2.7± 1.5 

4.3±0.7 

3.2± 0.6 

2.1±0.6 

5.6± 1.8 

54030 

2.7± 0.9 

4.0± 0.6 

2.2± 0.6 

1.2± 0.5 

0.6± 1.6 

54037 

3.6± 1.0 

3.7±0.6 

2.2± 0.5 

1.2± 0.5 

1.4± 1.4 

54040 

1.5± 1.5 

3.4± 1.0 

0.7± 0.8 

0.4± 0.6 

0.8± 1.7 

54049 

3.6± 2.2 

3.8± 1.1 

2.9± 0.7 

1.4± 0.6 

1.6± 1.9 

54053 

3.0± 1.1 

5.1±0.7 

3.4± 0.6 

1.5±0.6 

3.6± 1.7 

54056 

2.6± 1.3 

4.2± 0.6 

0.6± 0.6 

-0.8± 0.6 

-1.7± 1.7 

54059 

3.2± 1.3 

4.1 ±0.7 

1.9± 0.6 

1.3±0.6 

6.7± 1.9 

54062 

3.7± 1.2 

3.4± 0.6 

1.6± 0.6 

0.8± 0.5 

-0.6± 1.7 

54357 

3.5± 1.0 

3.6± 0.6 

1.8±0.5 

1.2± 0.5 

1.9± 1.4 

54359 

2.3± 1.4 

3.7±0.7 

1.6± 0.6 

0.7± 0.7 

0.2± 2.0 

54362 

3.3± 1.1 

3.3±0.6 

1.7±0.6 

1.7± 0.5 

-1.4± 1.5 

54376 

2.1± 1.2 

4.1 ±0.7 

1.4± 0.6 

0.7± 0.5 

2.3± 1.8 

54382 

1.8± 1.1 

2.3±0.6 

0.9± 0.6 

1.2± 0.6 

1.8± 1.8 

54385 

1.8± 1.1 

3.1±0.6 

I.l±0.5 

0.6± 0.5 

0.3± 1.4 

54411 

3.5± 1.0 

3.4± 0.6 

2.2± 0.5 

0.8± 0.5 

0.7± 1.4 

54413 

2.0± 1.4 

3.6± 0.7 

2.7± 0.6 

0.8± 0.6 

1.2± 1.9 

54415 

4.2± 1.1 

3.4± 0.9 

1.6± 0.5 

1.0± 0.5 

2.4± 1.6 

54422 

3.3± 1.3 

2.9± 0.6 

1.6± 0.6 

0.7± 0.6 

0.5± 1.6 


in the five optical bands (black empty stars, as in Fig.|^. The emis¬ 
sion in the u’g’r’i’ filters is consistent with a spectrum that rises 
steeply with frequency. On the other hand, the z’ band magnitudes 
are affected by the strong Ha emission, which, at the peak, con¬ 
tributes to about 25% of the observed flux in that band, and this 
contribution increases as the luminosity of the power-law contin¬ 
uum declines. To better characterize the spectral evolution of the 
flare, we first fitted the u’g’r’i’ data points with a simple power- 
law Li, oc i'~°‘ (or, equivalently. Lx oc and we plot 


the evolution of the spectral index a in panel (d) of Fig. [5] The 
slopes that we obtain are slightly steeper than the expected 
law of a geometrically thin, optica lly thick (un-truncated) accre¬ 
tion disc I Shakura & Sunva^l 19731) . but are co nsistent with previ¬ 
ous optical observat ions of well-sampled TDEs jGezari et alj2012l : 
iFloloien et alj|201^ . The overall ‘colour’ evolution of the flare, as 
diagnosed by the Stripe-82 u’g’r’i’ light-curves, is mild, and rem¬ 
iniscent of the almost constant-te mperature evolution seen in most 
well-sampled TDE light-curves l lGezari et al.ll201 iFIoloien et alj 
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Figure 5. Panel a\ the time evolution of the nuclear luminosity is shown with the coloured circles (where each colour identifies one epoch of the SDSS Stripe 
82 observations). The black circles are the total optical luminosity of SDSS J0159+0033, shown here for reference, together with the estimated (constant) 
luminosity of the host galaxy (red horizontal dot-dashed line). The empty triangles show the X-ray (2-10 keV) luminosity measured by XMM-Newton (first 
epoch) and Chandra (second epoch). The black solid line is the best fit power-law evolution of the flare decay, including all data points (with slope of 
— 1.59 ± 0.05). The dashed line, which is not a ht to the data, shows the long-term decline expected in most TDE models. The empty star shows the 

approximate location of the peak in our “fiducial” TDE model (see text for details), while the empty square at late times is the corresponding optical luminosity 
of the power-law component fitted to the BOSS 2010 spectrum, as described in section lT^ and plotted in Eig.|6]with a purple dashed line. Panel b and c show 
the evolution of the best fit photospheric radius (in units of the Schwarzschild radius for a 10® Mq black hole) and temperature (in units of 10^ K), obtained 
by fitting the optical photometric SED with a simple black-body spectrum (filled circles). Panel d shows the best-fit power-law slope a, where Lx oc 
In all panels, the black solid line connecting the points is drawn to guide the eye. 


|2014) . In fact, this behaviour is a challenge to the simplest model 
of a ‘bare’ viscous accretion d isc evolution for t he de cay phase 
of tidal disruption events, and iGuillochon et alj ( l2014h have ar¬ 
gued that an extended, large scale reprocessing layer should be 
present to produce the observed (almost) achromatic evolution of 
TDE flares. Such a reprocessor should be expected on the ba¬ 
sis of the results of hydrodynamical simulations, that show in¬ 
deed large amounts of debris present at distances ranging from the 
outer disc of the bound material that slowly returns to pericenter, 
ro « to the size of the expanding shell 


of unbound material expelled from the disrupted portions of the 
star, which moves away from the black hole (see section|4]below). 
As a simple phenomenological test of the above scenario, we have 
also fitted the same spectral evolution data shown in fig. with 
a single-temperature black-body spectrum, assuming, for simplic¬ 
ity, an emitting area Aph = 47ri?ph. Panel (c) of Fig.|^shows the 
evolution of the best-fit photospheric temperature (in units of 10^ 
K), while panel (b) shows the time evolution of the photospheric 
radius of the emitting surface, J?ph (in units of the Schwarzschild 
radius for a 10 ®Mq black hole). Indeed, photospheric temperature 
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Figure 6. Time evolution of the nuclear optical SED during the flare of 
SDSS J0159+0033. For each of the first 16 epochs of the cleaned Stripe- 
82 light-cui've, spanning the first 5 rest-frame years, shown in Fig. we 
plot the monochromatic luminosity, XL\, in units of 10'^^ ergs~^, in the 
u’g’r’i’z’ bands as filled coloured circles. The colour-code is identical to 
that of Fig. 1^ each colour corresponds to a different Stripe-82 epoch, as 
marked on the right of each SED. As in Fig.|^ the empty black stars indicate 
the baseline (constant) host galaxy SED, derived by taking the weighted av¬ 
erage of the photometric data-points of the last season of Stripe-82 observa¬ 
tions (last ten epochs. Autumn 2007) and subtracting the small power-law 
contribution that was estimated from the 2010 spectrum. The light grey line 
shows the difference spectrum between the peak flare (2010) and the host, 
with the corresponding difference-photometry data-points marked as blue 
solid circles (see section[^for more details). 


changes throughout the flare appear to be mild, if not negligible: 
the flare displays both temporal and spectral evolution very simi¬ 
lar to known UV/optical tidal disruption flares. We discuss in the 
following section the constraints on the TDE model parameters we 
can obtain by modelling these data. 


4 THE FLARE AS A TIDAL DISRUPTION EVENT 

A star of mass M* and radius R, that happens to reach a distance 
from a supermassive black hole of mass Mbh of the order of the 
star’s tidal radius Rt i?,will not survive the 
encounter unscathed. Depending on the star’s inner structure and 
on the penetration factor (3 = Rt /Rp (where Rp is the pericen- 
ter radius), various degrees of disruption will be unavoidable, with 
deep encounters (/3 > 1) causing the star to be comple tely torn 
apart jAval et alj200d : lGuillochon & Ramirez-Ruizl2013h . We can 
rewrite the tidal radius in units of the black hole’s Schwarzschild 


radius (Rs = 2GMBti/c^y 


Rt/Rs ^ 5.06 


f M,\ 


Mbh 
10’^ Mq 



( 1 ) 


If the black hole is massive enough, or the star compact enough, 
then we expect tidal disruption will only take place inside the event 
horizon (located at Rs for non-spinning black holes, and at f?s/2 
for maximally spinning Kerr black holes), and will be invisible to 
the outside o bservers. Mor e accurate General Relativistic calcula¬ 
tion (see e.g. lKesd^l2012h will then be needed to predict the inci¬ 
dence of TDEs from spinning black holes. 

In tidal disruption events, whatever fraction of the star’s mass 
is shed by black hole tidal forces, about half of it will remain bound 
to the central black hole. After an initial short ‘fall-back’ time, 
given by the time it takes the most bound material to re turn to peri- 
center assuming the s tar is initially on a parabolic orbit jRee si 19881: 
iMacLeod et al. f2012h : 


ffb ~ 0.37/3-^ 



Mbh 


1/2 



3/2 

yrs, 


( 2 ) 


the stellar debris will pile up near pericenter. A number of pos¬ 
sible mechanisms (general relativistic precession, hydrodynam- 
ical dissipation, comp ressive magneto-rotational instability, see 
iGuillochon et alj[2014l) could then lead to rapid dissipation of the 
material’s internal energy, circularizing the orbit of the debris. Vis¬ 
cous dissipation will then lead to accretion onto the black hole. 
Provided that both dissipation and viscous transport of angular mo¬ 
mentum can act efficiently, the rate of mass accretion onto the black 
hole is then fixed by the rate of mass return to peri center. This , 
in turn, can be computed using Kepler’s third law iRee iGUi), 
and depends critically on the specific binding energy of the stellar 
material at the time of disruption (dM/dE). For flat distributions 
{dMIdE = constant), this leads to the well-known evolution 
of the accretion rate at pericenter. Detailed numerical simulations 
of TDEs have mostly confirmed that this is indeed to be expected, at 
least in a ‘bolometri c’ sense, and at times longer than the peak ac¬ 
cretio n time (see e.g. Rannrez;Ruiz&Rosswog|200^ Lodatoe^^ 


I 2 OO 9 I : IGu illochon & Ramirez-RuizI 2013 : Guillochon et alj 2014 ). 

In fact. IGuillochon & Ramirez-RuizI 12013h have also performed 
a detailed parameter study of the expected light-curve of differ¬ 
ent TDEs as a function of black hole mass, stellar mass, penetra¬ 
tion factor and the star’s adiabatic index. In particular, they have 
shown how the peak accretion time and rate, as well as the asymp¬ 
totic slope, noo, of the accretion rate light-curve can all be used to 
put constraints on the parameters of t he encounter, including the 
masses of the star and the hole, and /3 IGuillochon et al.|[2Q14l) . In 
the following, we will use the scaling derived from these parametric 
study, but we stress that the issues of debris circularization, and the 
efficiency with which the rate of matter returning at pericenter ma¬ 
terial can be converted into an accretion rate onto th e black hole are 
currently open ones, and subject to a intense study iH avasaki et alj 

I 2 OI 5 I: IShiokawa et aLlfioT^: IGuillochon & Ramirez-Ruizll2015ah 


4.1 Constraining the parameters of TDE models 

Based on the outcome of hydrodynamical simul ation of (Newto¬ 
nian) TDEs, IGuillochon & Ramirez-Rui3 120131) have presented 
scaling relations for the most important light-curve parameters, as 
a function of black hole mass, stellar mass and radius, polytropic 
index of the stellar structure and penetration factor. In particular. 
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they showed that the peak accretion rate scales as: 


M, 


peak 


= aJ 


Mb 


V106 Mq 


- 1/2 


M, 


R 


(^) M,/y 


R. 


- 3/2 


while the peak time of the accretion rate follows: 

/ A/f \ / r? \ ®/2 




3i. 

Rq 


yrs. 


(3) 


(4) 


In the above expressions, A-, and represent rational func¬ 
tions of the penetration factor j3, evaluated from the numeri¬ 
cal simulations for different values of the polytropic index 7 
(assumed to be equal to 4/3 and 5/3 for high- and low-mass 
stars, respectively). Their form can be fo und in Eqs. A5-A8 of 
lOuillochon & Ramirez-Rui3 l l2013l[^15bh . To further simplify the 
analysis, we assu me a fixed mass-ra dius relation for main-sequence 
stars, as given bv lTout et alj l ll996l) : we fix the boundary between 
low- and high-mass stars at 0.6Mq, and the black hole mass to 
Mbh = 10® Mq, so that the overall family of TDE light-curve 
only depends on the penetration factor /3 and the star’s mass M*. 

The allowed region in this two-dimensional parameter space 
is determined by imposing two condition|3 the first is that the ac¬ 
cretion rate at the peak is sufficiently high to explain the observed 
peak luminosity of the nucleus, Topt.peak 1.1 X erg/s, (for 
a given radiative efficiency, e and optical-to-bolometric correction 

^opt)' 


Mpeak ^ 9.8 X 10-= (^) (^) ' Mo/yr. (5) 

The second is that the pericenter passage must occur outside the 
event horizon. For simplicity, we assume a non-spinning black 
hole, thus Rp/Rs = /3 ~^Rt/R s ^ 1 (see Eq.[T]above). For Kerr 
black holes, closer encounters are possible, and for maximally spin¬ 
ning black holes, the condition would read Rp > Rg = Rs /2. 

The left panel of Fig. Q shows the allowed range of the two- 
dimensional parameter space defined by the penetration factor /3 
and the star’s mass (in solar units), given the above constraints. 
We obtain a minimum mass of about = I. 2 M 0 and 

a quite narrow range of penetration factors, which widens up as 
more massive stars are considered. Just as a reference point for 
discussion, we mark with a black star in Fig. |7] a ‘fiducial’ pa¬ 
rameter combination allowed by the data: a 2 solar mass star 
with pericenter passage just 1.2 times smaller than its tidal ra¬ 
dius. Such an event would produce a flare of optical luminosity 
Topt.peak — 2.8 X 10'*"‘(e/0.1)(Kopt/5)-^, almost three times 
higher than the maximum caught by SDSS in 2000. The peak time 
would be «0.96 years after disrupti on, and, from Eqs. (A3-A4) of 
lOuillochon & Ramirez-Rui3 ( 1201 3l) . would evolve to the asymp¬ 
totic decay power-law slope Uoo ~ 1.64, consistent, within the 
uncertainties, with the observed value. Based on the same fidu¬ 
cial model, the amount of mass lost by the star would be AM « 
0.36Mq (« 18% of the star’s mass). Indeed, as a final consistency 
check, we derived the (bolometric) fluence of the flare, measured 
with only the available data points: J^boi « 4 x 10®=(Kopt/5) ergs. 
Assuming a 10% radiative efficiency, this corresponds to about 0.2 
Mq accreted onto the central super-massive black hole, consistent 
with the expectations of the fiducial TDE model described above. 


® Because the SDSS Stripe-82 light-curve did not well sample the rising 
phase of the outburst, and we do not know the exact time of the star’s pas¬ 
sage at pericenter, we cannot impose a third constraint on the peak time. 
As a reference, the time elapsed between the first and the second SDSS 
photometric data points is approximately 1.176 rest-frame years. 


The right panel of Fig. |7] shows the relation between the crit¬ 
ical black hole mass MsH.crit, above which stars are swallowed 
whole and no tidal disruption flare can be observed and the star’s 
mass, for the case of non-spinning black holes and maximally spin¬ 
ning black holes. 

The above estimates serve to emphasize, first of all, that a TDE 
explanation for the observed flare is energetically viable, and, sec¬ 
ondly, that well sampled light-curves of tidal disruptions flares by 
black holes of known mass could provide tight constraints on ac¬ 
cretion parameters that are typically elusive in steady accreting sys¬ 
tems, such as the radiative efficiency (and thus th e black hole spin), 
and t he accretion flow bolometric corrections dGuillochon et alj 
l2014h . 

To conclude this section, we note that a prominent re¬ 
brightening flare was observed about 4 rest-frame years after the 
peak. This sub-flare is characterized by a relatively constant con¬ 
tinuum slope and photospheric temperature, and has a fluence 
of about 5.5 x 10®^(Kopt/5) ergs, corresponding to 14% of the 
fluence of the whole flare. Among well-m onitored TDE so far , 
only the 5w//i-selected ‘relativistic’ ev ents dBurrows et al. IH; 
iBloom et al.ll201 ll : ICenko et alj|2012bh did show clear structures 
in their declining light-curves, but it is possible that this is simply a 
selection effect due to lack of long-term dense monitoring of more 
recently discovered optical TDEs. Lacking a large sample of well 
sampled optical light-curves of TDEs, it is beyond the scope of this 
paper to speculate about the nature and causes of this apparent sec¬ 
ondary flare event, but one should consider testing hydrodynamical 
models of tidal disruptions, inc luding the expectations from binary 
events dMandel & Levin|[2015h . to assess whether such events may 
pose a real challenge to the TDE interpretation of this flare, and 
point towards an alternative explanation in terms of AGN accretion 
physics, as we discuss below. 


5 DISCUSSION 

A pre vious search for TDEs in the Stripe-82 data dvan Velzen et alj 
1 201 ih did not select SDSS J0159-I-0033 as a good tidal disruption 
candidate. We believe there are two concurrent reasons for this: first 
of all, the flare we observe is rat her slow (taking more tha n 5 years 
to return to quiescence), whilst Ivan Velzen et alj ( 1201 ih selected 
preferentially flares which lasted only one season o f Stripe-82 ob¬ 
servations. More critically. Ivan Velzen et ^ l20ll|) excluded from 
their final list of TDE candidates all those with known AGN spec¬ 
troscopic classification. It is just due to a fortuitous coincidence 
that the first 2000 SDSS spectrum, showing the prominent broad 
emission lines that led to type-1 AGN classification, was taken very 
close in time to the flare peak, and we suspect this crucial detail had 
been overlooked. 

Of course, the main question we would like to answer is the 
following: is the observed flare from SDSS J0159-(-0033 the result 
of the tidal disruption of a star, or just a particularly strong and rapid 
AGN outburst? To a certain extent, this is an ill-posed question, as 
TDE are in fact just a specific class of AGN outbursts. The anal¬ 
ysis presented in the previous section, however, suggests that the 
observed light-curve evolution, and, in particular, the overall ener¬ 
getics of the event are consistent with the amount of accreted mass 
being a substantial fraction of a massive main-sequence star. In gen¬ 
eral terms, well-sampled light-curves are probably the most pow¬ 
erful tool in order to distinguish TDE-induced flares from stochas¬ 
tic AGN variability, as demonstrated by the recent example of the 
nearby source IC 3599 jCampana et alfcoT^ . 
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Figure 7. Left: The allowed range of the TDE model parameters (3 = (penetration factor) and M* (star’s mass) is shown as a doubly shaded (blue and 

red) wedge. It is derived by the combination of two independent constraints: (i) pericenter passage must occur outside the SMBH event horizon (i?p > i?g, 
thick blue line, for the two extreme cases of non-spinning black hole and thin blue line for the case of maximally spinning hole, Rp > -Rg); and (ii) the 
peak luminosity of the flare has to be at least as high as the observed one (assuming a optical-to-bolometric correction of 5 and a radiative efficiency of 10%, 
solid red line). The black dashed lines show the contours of equal peak time, in rest-frame years. The black star mai'k a ‘fiducial’ model, discussed in the 
text, for illustration puiposes. Right: The relation between the critical black hole mass MBH,crit (x-axis), above which stars are swallowed whole and no 
tidal disruption flare can be obseiwed and the stain’s mass (y-axis) is shown as a black solid line for the case of non-spinning black holes (leftmost line) and 
maximally spinning black holes (rightmost line). The vertical dashed area mark the range of possible SMBH masses estimated for SDSS J0159+0033 on the 
basis of the observed broad emission lines. 


Irrespective of the source of the accreted material (star or gas 
cloud, for example), it is clear that it must have been deposited very 
close to the central black hole. If the accretion energy we se e was 
generated in a standard, viscous. [shakura & Sunvae^ jl973h disc, 
the viscous (or infall) time, can be written as: 

tvisc 1.7 X 10“'‘m8ao.i ® 

where mg is the BH mass in units of 10 ®Mq, ao.i is the vis¬ 
cosity parameter in units of 0.1 and r is the radius of the disc in 
units of Rg. The disc scale-height in the radiation-pressure dom¬ 
inated inner part of the disc, which is the region relevant for the 
luminous accretion episode discussed here, can be expressed as 
H/R ~ 20\EddJ{r)/r, where AEdd = T-boi/T/Edd is the Ed¬ 
dington ratio (AEdd ~ 0.04(«:opt/5) at the observed peak, and a 
factor of three higher at the fiducial peak) and J(r) = 1 — yrn^/r 
is the Newtonian correction factor from the no-torque inner bound¬ 
ary condition. Thus, viscous times shorter than one year require 
disc sizes smaller than 10 gravitational radii (i.e. smaller than a few 
times the pericenter distance, with Rp « 3.37?g in our fiducial 
model), with the exact details dependent on the disc inner bound¬ 
ary condition (and BH spin). Progressive draining of a thin disc 
substantially larger than this is plainly incompatible with the ob¬ 
servations. However, thermal instabilities can also operate in the 
inner parts of the disc, which will be faster by a factor of the or¬ 
der (R/H)^. In fact, there are a number of potential mechanisms 
that theories of accretion discs have put forward to explain very 
rapid (faster than viscous) variability in the inner region of an ac¬ 
cretion discs, such as thermal instabilities, large scale waves in the 


inner accret ion disk, hydromagnetic winds, reproce ssing of UV or 
X-ray light l lKrolik et alj2005l : TLaMassa et alj20l3 . and references 
therein). 

The X-ray spectrum observed by XMM-Newton close to the 
peak of the lightcurve is consistent with a power-law of index 2.1 
jLaMassa et alj|201^ . similar to what typically observed in lumi¬ 
nous AGN. Simple thermal models for the emission from a TDE do 
not predict such a ’corona-like’ spectrum, very much as the stan¬ 
dard Shakura-Sunyaev theory of geometrically thin and optically 
thick accretion discs does not predict the ubiquitous Comptonis- 
ing medium inferred from X-ray spectra of AGN. Indeed, most 
TDE candidates detected by ROSAT, XMM-Newton and/or GALEX 
showed mostly very soft spectra, well-fit by a black b ody model 


and/or much steeper nower-law (Brandt et al 

I 995 I; Bade et al 

19961; 

Komossa & Greineilll999bl:lGreiner et al. 

2000l;lGezari et al 

2008|; 

Esauei et alj 2008h. But harder X-ray spectra have been ob- 


served recently from candidate TDEs, too dCenko et al.l 12012al : 
iNikolaiuk & Waited l2013l : ISaxton et al.ll2014h . In general, too lit¬ 
tle is currently known about the spectral formation mechanisms in 
TDE, across the entire electromagnetic spectrum, to be able to re¬ 
ject the TDE interpretation on the basis of the X-ray spectral shape 
only. Moreover, we note that, in general, optically selected TDE 
tend to have low photospheric temperatures and be X-ray faint, 
while X-ray selected TDE tend ed to be hot (when thermal) and 
optically faint (see e.g. Eig. 4 in lGezarill2012h . clearly suggesting 
we are barely scratching the surface in our understanding of TDE 
selection effects. 

There are, moreover, other pieces of evidence that are not 
straightforward to interpret within the simplest TDE scenario we 
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have described so far. In particular, the properties of the emission 
lines observed in the spectra are not easily explained by a TDE 
model, as we discuss in the following subsections. 


5.1 The broad line region 

Only since the discov e ry of the first op tically selected TDEs 
dvan Velzen et al. 1 1 201 ll : iGezari et alj 1201^ . and the consequent 
availability of early spectroscopic follow-up observations, it has 
become evident that the emission from tidal disruptions of stars 
is accompanied by the presence of broad emission lines in their op¬ 
tical spectra. The broad emission line phenomenology is diverse, 
with some objects showing only high-ionization lines (typically 
Hell A4 686), others only B almer lines, and others still both He and 
H lines dArcavi et al.ll2014h . Measured line widths vary between a 
few thousands to about ten thousand kms“^. If these broad emis¬ 
sion lines are produced by stellar debris illuminated by the tidal 
disruption flare, we can immediately rule out them being located at 
a distance from the hole of the order of the circularization radius 
(typically assumed to be equal to 2Rp), as this would imply (Kep- 
lerian) velocity widths of a few times km s“^. In fact, two pos¬ 
sible alternative si tes of broad emission lin e production have been 
invoked for TDE. iGuillochon et alj d2Q14t) argues that the bound 
material extending out to ro (see section [ATI above! will constitute 
an elliptical accretion disc, growing inside-out, eventually reaching 
densities and ionization states leading to the generation of broad, 
permitted atomic emission lines, very much like those generated 
by long-lived AGN. According to this scenario, high-ionization 
lines should appear first, followed by lower ionization lin e s later 
on during the outburst. Alternatively, IStrubbe & Ouataerl ( l2009l) 
considered the outbound stellar material, flying away from the dis¬ 
ruption site and extending out to i?max = ~ 

2.2 X — ip)/2/f] cm, as the site of line emis¬ 

sion processes. They, however, predicted both emission and ab¬ 
sorption lines, and significant bulk blue-shift of the lines, con- 
tr ary to what is typically ob served. Moreover, as discussed in detail 
in IGuillochon et al] ( l2014ll . Newtonian hydrodynamic simulations 
clearly show that the unbound material will be confined to a thin 
filam ent, covering a sm all solid angle as seen by the central black 
hole jKochanekl 1 1 994 ) . hardly large enough to produce emission 
lines with the observed equivalent width. 

The broad Balmer emission lines (Ha, H/3 and H 7 ) observed 
in the 2000 SDSS spectrum of SDSS J0159-I-0033 do appear at 
first sight like ordinary QSO broad lines: they are not strongly 
asymmetric, have no absorption through or P-Cygni profile, and 
have measured FWHM ranging between 3.4 and 4.5 x 10® km s~® 
(see Table [T]l. If they are produced by a mechanism analogous to 
that which generates the broad line region (BLR) in most AGN 
(an assumption we have implicitly used by adopting the “single 
epoch” virial black hole mass estimates for this object), then they 
should follow the BLR luminosity-size relation (7?blr oc L^, with 
S = 0.4 — 0.6), which woul d put them at a distance of about 30 
light days ( iBentz et aD l2006h . i.e. « 2.6 x 10 ® fls for a 10 ® 
centra l black holo This is inconsistent with the IGuillochon et ^ 
( I 2 OI 4 I) ‘ inner accretion disc’ interpretation, as ro is only about half 
this value at the time of the 2000 SDSS spectroscopic observation, 


® Similar values for the BLR size would be obtained by adopti ng the empir¬ 
ical sc aling between X-ray 2-10 keV luminosity and BLR size. lGreene et alj 

J2010l) . 
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and only reaches these distances about three rest-frame years after 
the observed peak. 

However, the most serious problem for the interpretation of 
the broad emission lines as originating in the stellar debris, are 
the implied gas mass by the standard ionization/recombination 
model for BLR. Comparing the measured equivalent widths of 
Ha and H/3 in the 200 0 SDSS spectrum with those predicted by 
iKorista & Goad l l2004l) as a function of the ionizing flux and BLR 
particle density, we derive a particle density of about log n « 10.5 
(assuming an ionizing luminosity of 10^^ergs“®, and a distance 
of the BLR of 30 light days). The amount of material needed, 
assuming spherical geometry is a few hundred solar masses (see 
iBaldwin et al.ll20()3h . obviously much more than could have pro¬ 
vided by the disrupted star. 

For the particle densities inferred for the BLR emitter, re¬ 
combination is very fast, so the lines respond effectively instan¬ 
taneously to flux variations from the central source. Thus, the rapid 
luminous flare we observed between 2000 and 2005, must have il¬ 
luminated a pre-existing structure, whose spatial distribution and 
kinematics closely resemble those of BLR of actively growing 
SMBH. Further evidence for this can be derived by the fact, no¬ 
ticed bv iLaMassa et alj ( l2015h . that the measured Ha FWHM be¬ 
tween 2000 and 2010 scales as the continuum luminosity to the 
— 1/4 power, as expected if gas in the BLR moves with Keplerian 
velocity. 

The likelihood of the observed configuration, with the broad 
emission lines appearing transiently in response to a rapid nuclear 
TDE, depends on the lifetime of the broad-line region structure in 
objects where the central black has ceased accreting at substantial 
rates. This must depend, in turn, on the dynamical status of the 
BLR itself (is it outflowing, inflowing, circulating?), and, critically, 
on the time elapsed since the last major accretion episode onto the 
SMBH. A detailed analysis of these issues is beyond the scope of 
this paper. In any case, the conclusion that such a configura tion is 
possible was reached recently also by IPennev et alj ( l2014ll . who 
studied the peculiar, decades-long flare of the central black hole 
in Mrk590, which showed dramatic brightening and subsequent 
disappearance of strong broad emission lines over little more than 
20 years. 

5.2 The narrow line region 

It is not just the properties of the observed broad emission lines 
that suggest the flare we observed in SDSS J0159-I-0033 occurred 
in a formerly active galaxy. The system in question is massive and 
strongly star-forming, as indicated both by the SED fitting and by 
the strong [O II] emission line detected (see section 12.21 and Ta¬ 
ble [T] above). In the n earby universe, these ar e the systems most 
likely to host an AGN jKauffmann et al.ll2003l) . The narrow emis- 
sion line diagnostics (see Fig. O reveal that the emitting material 
was probably ionized by a source harder than that associated with 
pure star-formation. AGN activity would of course provide a suffi¬ 
ciently hard ionizing source, but the exact amount of any putative 
AGN contribution to the observed narrow line emission is harder to 
disentangle. SDSS J0159-I-0033 appears to be located close to the 
boundary of the so-called ‘transition’ region of parameter space, 
between the objects whose line ratios imply stellar sources domi¬ 
nate the ionization, and those for w hich AGN emission is n eeded 
to produce the observed line ratios. iKauffmann et al.l l l2003l) esti¬ 
mate that the AGN fractional contribution for objects in the transi¬ 
tion region ranges between 30 and 90 percent. In the most extreme 
case, then we can estimate the maximum luminosity that any pu- 
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tative AGN must have had in the past to produce 90 percent of the 
observed [Olll] flux (corresponding to -£/[oiii] — 3 — 5 x 10^^ 
ergs"i). 

There is a large body of literature on the correlation be¬ 
tween [O III] and hard X-ray luminosity in AGN, with somewhat 
contrasting results on the exact numerical value of the correla¬ 
tion index and on the mean luminosity ratio. In general the ra¬ 
tio of observed (not extinction corrected) [O III] to X-ray (2-10 
keV) luminosity ranges between 0.02 and 0.08, depending on the 
sample selection and source classificatio n i Mulchae^_et_^ 199^ 
Heckman et all 1 20051 : iTrouille & Bargej 120101 : de Gasperin et all 


201 ih . Thus, we can conservatively say that the central AGN must 
have had an X-ray luminosity L 2 —lokev s; 2 X ergs"i (i.e. 
equal or smaller than the X-ray luminosity measured by Chandra 
in 2005, Z/2-iokeV = 1.2 x 10“^^ erg s“^) about one narrow-line re¬ 
gion light-crossing time ago. Us ing the empirical relation between 
NLR size and [O III] luminosity jGreene et alj20li] : lHainline et alj 
l2013h . we derive a NLR size of about 4 kpc, corresponding to a 
look-back time of about 13 thousand years. As discussed in sec¬ 
tion |%T] the historical optical data pre-flare (i.e. before 1998) are 
consistent with no strong AGN emission being present above the 
host galaxy stellar emission, but the uncertainty on the measure¬ 
ments is such that a low-level nuclear activity cannot be ruled out 
either, so that it is not possible to firmly constrain the possible AGN 
activity level in the galaxy prior to flare discussed here. 

Powerful TDE, just like longer episodes of nuclear activ¬ 
ity powered by gas accreti on, should imprin t in the host ISM an 
ionization ‘echo’ (see e.g. lYang et ^l20I3tl detectable as high- 
ionization narrow emission lines, variable on timescales longer than 
the light-crossing time at the location of their production. Variabil¬ 
ity of strong coronal emission lines, possibly associated to the fad- 
ing echo of a nuclear tidal disrupt i on ev e nt has been repor ted by 
iKomossa et aP l l2009tl : lmng et alj ( l2Ql3l : lYang et al] ( l2013ll . Sys¬ 
tematic spatially-resolved spectroscopic surveys with integral field 
unit of nearby galaxies (such as SAMI or MaNGA. ICroom et alj 
l20I2l:lBundvetalj20I5l) could be able to detect signatures of such 
echos over the entire body of nearby galaxies, thus probing up to 
timescales of the order fecho ~ Rbb/c ~ 3 x 10^(i?eff/10kpc) 
yrs. Given current estimates of the intrinsic rate of tidal disruption 
events per galaxy Ttdb = l/lrec ~ 10“® yr~^ (where we have 
introduced the recurrence time free; see sections [T] above and l53] 
below, and references therein), it is possible that a large fraction 
of nearby galaxies, of the order of lecho/lrec ~ 30%, could be 
found displaying off-center high-ionization emission line regions, 
thus probing the triggering properties of nuclear black holes, and of 
tidal disruption events (which have the shortest recurrence times) in 
particular. 


5.3 Tidal disruption flares vs. gas accretion in galactic nuclei 

We discuss here a simple unified model for the stochastic ‘acti¬ 
vation’ of a nuclear black hole, be it via tidal disruption flares 
(‘TDE’) or accretion of gas parcels from the galactic ISM (‘ACC’). 
Without loss of generality, and for the sake of simplicity, we as¬ 
sume here that individual accretion events are all characterized 
by a light-curve composed of two main phases: a ‘peak’ phase, 
where the luminosity rises rapidly and saturates at a peak value 
for a time TijMsH, z) (where, here and in the following, the index 
i = ACC, TDE, depending on the case we are interested in), and a 
‘decline’ phase, where the light-curve follows a power-law decline 
L oc t ’’L While this is now a well tested theoretical framework 
for TDE, as we have discussed above, it has also been shown by a 


number of authors that AGN light-curves of this sort do indeed re¬ 
produce generic features of the active galactic nuclei distributions 
such as luminosity functions, accretion rate distributions, cluster¬ 
ing properties, etc. 1 Yu & Lull2004l : iHopkins et alj|200^ : IYu & Lul 
iTOoklSonoli et al.l2009l) . We further define the rates of nuclear ac¬ 
tivation as ri(MBH, z). It is straightforward to show that, under 
these simplified circumstances, the probability of a galaxy hosting 
a black hole of mass Mbh which has bolometric luminosity be¬ 
tween Lboi and Lboi + d log Lboi is given by: 


Pi{Lhol) — 

li 




-l/7i 


\ -^'peak J 


/(Ebol/Epeak), (7) 


where /(Lboi/ipeak) is a ’cut-off’ function which is « unity for 
Lboi Tipeak, and rapidly declines for Lboi > Lpeak (typically 
exponentially). 

In the case of tidal disruption of stars {i = TDE), we can 
write an approximate solution neglecting any redshift dependence, 
fixing 7tde = 5/3, ttdb = fpeak( Mbh, M*,/3) (Eq. |4j, and 
Lpeak = £Mpeak(MBH,M*,/3)c^ (Eq. |3j. Eor simplicity, we 
discuss here only tidal disruptions of main-sequence stars, and 
ref er the reader to the m ore accurate, and comprehensive work 
of iMacLeod et alj ( l2013h . who have computed the contribution 
to PTDB(Lboi) of giant stars, either disrupted, or undergoing re¬ 
peated, episodic, mass-transfer events (‘spoon feeding’), which 
tend to dominate at very low luminosities. As we discussed in the 
introduction, there is substantial uncertainty in the true rate of tidal 
stellar disruption in galactic nuclei, with current observationally- 
determined rates lying factors of a few to ten below the theo¬ 
retical expectations. Here, for illustration purposes, we adopt a 
generic (i.e. independent of black hole mass) rate of Ttdb = 
3 X but consider, in the following, a one-order- 

of-magnitude allowed range around such a value. The resulting 
probability distributions, Ptdb (Lboi), are shown with a solid black 
line in Eig.[^ for three different values of Mbh. 

In the case of gaseous accretion, most attempts to reproduce 
the AGN luminosity function have implied peak luminosities of 
the order of the Eddington luminosity. However, both the trigger¬ 
ing rate rAcc(MBH, z) and the peak duration tacc(Mbh, z) are 
essentially unknown. In fact, it is unlikely that a single mech¬ 
anism could provide a complete description of (gaseous) AGN 
triggering. Nevertheless, multi-wavelength surveys of AGN have 
been used, in recent years, to provide empirical constraints on 
PACc(Lboi), by studying in a uniform way the properties of the 
host galaxies of complete AGN samples typicall y selected in the 
X-ray band, to minimi ze the ‘obscuration’ bias |Bongionioet_^ 
12012 : lAird et al] l2012h . Interestingly enough, lAird et all (1201 2h 


showed that the PACc(Lboi) is independent of black hole mass, 
but strongly evolving with redshift p (Lboi) oc (1 -I- z)^'^, a result 
confirmed bv iBongiomo et all (l2012h . both qualitatively and quan¬ 
titative ly. We make use here of the analytic solution of lAird et alj 
( l2012h . who expressed PACc(Lboi) as a broken power-law, with 
slope I/ 7 ACC = 0.65 and a rapid decline above the Eddington 
luminosity. 

The very fact that, on average, the overall AGN population is 
consistent with 7acc — 1.5 is interesting, as this is not only a 
steeper decline than predicted by theoretica l models of viscously 
evolving accret ion discs (jacc ~ IT — 1-2, ICannizzo et al.ll9^ : 
IYu & Lj|2004h . but also shallower than the ‘self-regulated’ light- 
cu rve of the quas a r-feed back dominated major mergers simulated 
by iHopkins et alj (l2005l). which pred ict 7acc ~ 2. Instead, the 
empirical results of Aird et all l l2012ll indicate that, in an average 
sense, the light-curves of AGN do resemble those of tidal disruption 
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Figure 8. The probability for a galaxy hosting a black hole of mass Mbh to have bolometric luminosity in the range Lboi to Lbol + d log Lboi are shown for 
three values of Mbh (10®Mq, left; lO^ Mf?), center; 1Q® M(T), right). In each panel the three red lines represent the overall AGN probability function derived 
empirically from AGN X-ray surveys bv lAird et al.1 20lj), f or thre e redshift values {z = 0, solid; z = 0.3, dashed; and 2 : = 0.6, dotted), the black solid line 
is the rate for TDE computed averaging Eq.[3over a iKroupal bOOlh IMF from 0.3 to 30 Mq . The grey band represents a one-order-of-magnitude uncertainty 
in the TDE rates. In each plot the thin blue solid line shows the ratio of the black to red solid lines, i.e. the z = 0 fraction of luminous galactic nuclei powered 
by tidal disruptions of main-sequence stars. 


flares (we have 7tdb — 1-66 and 7acc — 1-5); the difference 
in the overall probability functions between the two populations 
must then be the outcome of widely different triggering rates, event 
durations and (to a lesser extent) peak luminosities. 

Fig. [8] shows, for three different values of Mbh the prob¬ 
ability of hosting an active nucleus with luminosity Lboi (per 
unit logarithmic lu minosity int e rval) f or the full AGN population, 
as parametrized by lAird et all l l2012h (red lines), and for main- 
sequence stars TDE (solid black line with grey b and uncertainty), 
where we have averaged over a iKroupal <200ih IMF, and taken 
/3 = min[l, ,gniax(MBH, M*)], where ,5max(MBH, M») is the 
maximum value of the penetration factor still compatible with a 
tidal disruption (i.e. that value for which Rt/Rs ~ 1, see Eq. [T] 
above). The plot also shows, in each panel, the ratio of the TDE 
probability function to the z = 0 AGN probability function, which 
represent the fraction of active nuclei (of a given mass) at any time 
which are in fact powered by the disruption of a main sequence star. 

We can also weigh the probability functions Pi{L) with the 
local black hole mass function (BHMF) and compute the total frac¬ 
tion of TDE-powered nuclei as a function of nuclear bolometric lu¬ 
minosity for a volume limited sam ple of local AGN. We adopt the 
BHMF of iMerloni & Heinj ( l2008h . and show the results in Fig. 

A similar calculation was performed bv iMilosavlievic et al.l l l2006r) . 
who compared TDE volumetric rates with X-ray AGN luminosity 
functions, and reached similar conclusions as those we can draw 
here: the fraction of TDE in current AGN samples is highly uncer¬ 
tain, due to the still poorly known intrinsic rates, but, in the range 
10'*^ < Lboi < 10^"^ergs“^ could be as high as 10%. The com¬ 
bined effect of the Eddington limit and of the fact that bigger black 
holes would swallow most main-sequence stars whole, produces a 
rapid drop of the TDE fraction at bolometric luminosities above 


lO'^^ergs”^, i.e. into the quasar regime. The exact behaviour of 
the TDE fraction at such high luminosities, however, critically de¬ 
pends on the SMBH spin distribution, as rapidly spinning black 
holes can tidally disrupt stars even if they are very massive (see the 
right panel of Fig.[7]above). On the other hand, a proper general- 
relativistic treatment of TDE could modify the p(Lboi) even at low 
Mbh, where the Lense-Thirring precession of the stellar debris 
may cause severe de lays in the onset of accretion, increa sing fpeak 
and reducing Lpeak ( iGuillochon & Ramirez-Rui3l2015all . 


As a term of reference, we also show in Fig. [9] the observed 
fraction of all X-ray selected AGN in a highly spectroscopically 
complete sample in the redshift range 0.1 < z < 0.6 (from the 
BOSS follow-up of the XMM-XXL survey, Menzel et ah, in prep.), 
which are optically classified as ’elusive’, or ’X-ray Bright and 
Optically Normal Galaxies’ (XBONGS), i.e . lacking any spectro¬ 
scopic signature of AGN emission (see e .g. IComastri et alj|2002 : 
ISmith et S] 12014 IPons & WatsonI I20l4) . At these luminosities, 
about 6-20% of all X-ray selected AGN are ’elusive’, and we sug¬ 
gest that a non-negligible fraction of them are powered by TDE 
which have faded away by the time the optically spectroscopic 
follow-up was performed (3 to 10 years after the X-ray data were 
taken in this field). Further clues on the nature of elusive AGN, 
and on the true fraction of TDEs among them could come from a 
detailed study of the IR emission from those sources, which could 
reveal either the presence of a deeply obscured AGN or the tran¬ 
sient echo of the TDE emission. Unfortunately, no systematic, deep 
NIR/MIR data are available for the full X-ray selected XMM-XXL 
sample. 
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Figure 9. The overall fraction of active galactic nuclei powered by TDE at 
z = 0 is shown as function of nuclear bolometric luminosity (black solid 
line), with a grey band representing a one-order-of-magnitude uncertainty 
in the TDE rates. The dashed and dotted lines show the evolution of the TDE 
fraction at 2 = 0.3 and 2 = 0.6, respectively, if the TDE rates are assumed 
to be unchang ed, while the tota l AGN triggering rate evolves as observed 
empirically bv lAird et al.N2012h (PACC oc (1 + 2 )® ®). The empty squares 
are the observed fractions of “elusive” AGN among all X-ray selected AGN 
in a highly spectroscopically complete sample (with flux Fg-iokeV > 2 X 
10“^^ ergcm~^ s“^) in the redshift range 0.1 < 2 < 0.6 from the BOSS 
follow-up of the XMM-XXL survey (Menzel et al., in prep.). 



6 CONCLUSIONS 

We have analysed in detail the light-curve and spectra of the flare 
observed between 1998 and 2005 in SDSS J0159-I-0033, a mas¬ 
sive star-forming galaxy in the SDSS Stripe 82 field. The most 
striking features of the system are: (i) a rapid increase in the X- 
ray and optical/UV luminosity of the source by at least one order 
of magnitude, and (ii) the dramatic change of its optical spectrum 
over a ten year timescale. The object changed from a typical type- 
1 AGN emission, with prominent broad emission lines of Ha and 
H/3 and a blue continuum over-imposed on a galaxy spectrum, to 
a galaxy-dominated continuum. This latter spectrum shows strong 
narrow emission lines, which could have been ionized by either 
newly formed stars or (more likely) by past AGN activity, and, in 
addition, a weak broad Ha emission line, suggesting low-level of 
ongoing activity a few years after the peak. 

We have studied the light-curve of the object over the last ~30 
years, in optical/UV and X-rays, including a well-sampled period 
coincident with the flare decay from the SDSS Stripe 82 survey. 
The variable nuclear emission, both in X-rays and in the optical 
bands, is well described by a single flaring episode, with fast rise 
and a power-law decay with exponent rioo = —1.59 ± 0.05, con¬ 
sistent with detai led predictions of hydro dynamical models of tidal 
disruption flares l lGuillochon et alj2014t) . The overall fluence of the 
event indeed suggests that the amount of accreted mass is a sub¬ 
stantial fraction (at least about 20%) of a massive star. The decay 
light-curve also shows clear structure, with a re-brightening flare 


observed about four rest-frame years after the peak, which has a 
fluence of about 14% of the total fluence of the event. 

We have used late-time photometric and spectroscopic obser¬ 
vations of the source to build a host galaxy model, that we used to 
carefully analyse the nuclear (host-subtracted) light-curve of the 
flare in all five SDSS bands. This has allowed us to study the 
SED evolution of the flare, which we found, consistent with most 
optically-selected TDEs observed to date, to display very little tem¬ 
perature (colour) evolution, as predicted by models in which the ac¬ 
cretion energy releas ed in the event is reproce ssed by dense, large- 
scaleheight material jGuillochon et al.ll2014h . Assuming standard 
(and conservative) bolometric corrections, this is among the most 
luminous non-beamed tidal disruption flares discovered so far, and 
the only one produced by a black hole as massive as 10® solar 
masses. Perhaps due to the serendipitous nature of its discovery, 
it is also one of the few cases in which both optical and X-ray 
emission contribute substantially to the bolometric luminosity of 
the event, suggesting that accretion onto black holes from tidally 
disrupted stellar debris can also give rise to powerful Comptonis- 
ing coronae. In most previous cases, optically selected TDE had 
low photospheric temperatures and were X-ray faint, while X-ray 
se lected TDE te nded to be hot and optically faint (see e.g. Fig. 4 
in iGezar il l2012l) . clearly suggesting we are barely scratching the 
surface in our understanding of TDE selection effects. 

If the continuum evolution and the overall energetics of the 
event strongly suggest a TDE interpretation, the emission line spec¬ 
tra tell a more complicated story. The broad line region, which re¬ 
sponded almost instantaneously to the strong illumination from the 
central flaring source in 2000, is too extended (about 30 light days) 
and massive (more than hundred solar masses) to have been cre¬ 
ated by the debris of a disrupted star, and must have been left, in 
the dark, from previous episodes of more prolonged accretion onto 
the central black hole. Also, the narrow emission lines, which we 
see unchanged between 2000 and 2010, possibly indicate that the 
central black hole was active as little as 10^' years ago, although 
it is not straightforward to estimate exactly at which level, as the 
narrow emission lines could be also partly excited by the powerful 
star-formation observed in this galaxy. 

Such a complex interplay must not be too exceptional, also 
given the fact that this TDE was discovered in a relatively small 
parent sample. Within a simple, common framework in which both 
TDE and gas accretion onto a SMBH are modelled as stochas¬ 
tic events with homologous light-curves, but widely differing trig¬ 
gering rates and durations, we showed how the available empir¬ 
ical and theoretical constraints suggest that TDE may represent 
a substantial fraction (between « 1% and 10%) of all putative 
AGN in “snapshot” surveys, at least for bolometric luminosities 
Lhoi < 10^® ergs“^. In particular, we expect that tidal disruption 
flares should be responsible for at least a fraction of those X-ray 
selected AGN, which, upon later follow-up analysis, do not reveal 
any evidence for AGN-induced signatures in their long-wavelength 
spectra (suc h as the so-called XBONGS, IComastri et al.l 1 20021 : 
ISmith et alj[2014h . Systematic re-observations of previously clas¬ 
sified broad l ine AGN has only been done for small samples of 
nearby AGN: [^ott et alj ( 120141) have indeed reported a consider¬ 
able n umber (8/97) of ‘changing look’ sources (see also lBarth et alj 
I 2 OI 5 L for the results of a recent AGN spectroscopic monitoring 
campaign). We are, indeed, used to think of AGN as being either 
type 1 (un-obscured) or type 2 (obscured), based on the properties 
of their optical emission lines. However, on timescales longer than 
the duration, r^, of episodes of high-luminosity (‘flares’), optical 
emission line classification is a transient variable. This is particu- 
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larly important for tidal disruption events, of course, which have 
flare durations of a few years, at most, but could be relevant for 
other forms of black hole fuelling, too. 

A few lessons can be learned from this work: probably the 
most compelling is that an accurate observational determination 
of TDE rates requires searches that do not a priori exclude ob¬ 
jects previously classified as AGN, and needs to fully account for 
the exact time when such classification was firstly made. Also, we 
showed that TDE confirmation is far from straightforward, and re¬ 
quires well sampled optical light-curves, and would be aided by 
repea ted, early spectroscopic observations, too. LSST dlvezic et alj 
l2008h . in combination with flexible, rapid responding spectroscopic 
surve y instruments (such as 4MOST, for example, Ide Jong et alj 
l2014h . will be instrumental for our understanding of nuclear vari¬ 
ability and AGN triggering. In addition, deep and wide X-ray ob¬ 
servation with the right cadence have the advantage that the TDE- 
produced light becomes confused with the host galaxy stellar back¬ 
ground at much later stages at these energies, thus allowing flares 
to be identified with less well-sampled light-curves and on longer 
timescales. Thanks to this fact, the upcoming eROSITA on SRG 
jMerloni et alj 2()12[) could yield as many a s two TDE candidates 
alerts per day l lKhabibullin & Sazonovll20l^ . requiring substantial 
dedicated resources to be mobilized in order to follow most of them 
up. 
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APPENDIX A: PHOTOMETRIC CALIBRATION OE DSS 
PALOMAR AND UK-SCHMIDT PHOTOGRAPHIC DATA 
AGAINST SDSS 

In order to convert and calibrate the pre-1998 photographic imaging data- 
points from the D SS (as shown in Fig. [T] from GSC 2.3.2 catalogue, 
iLasker et al.l2008ll into the SDSS photometric system, we cross-coiTelated 
point-like detections from the GSC 2.3.2 catalogue lying within 2 degrees 
from SDSS J0159+0033 with SDSS photometric objects (treated as a truth 
sample) in order to measure the bandpass conversion formulae and residual 
scatter for the photographic plate photometry. To relate the SDSS g’, r\ V 
bands to the observed j, F, N and V magnitudes, we used fitting formulae 
that are linear in SDSS colour and quadratic in the photographic magnitude, 
but the quadratic term turns out to be negligible in all but the F-band con¬ 
version. We obtained the following expressions for the SDSS magnitudes: 

g'ihg' - r') = 0.38811 + 0.96837j - 0.089848(g' - r' - 1) (Al) 

r'{F,g'-r') = -8.8646-l-1.4706F+0.10263(F-20)2-l-0.10345(3'-r' 

(A2) 

i'{N, r' - i') = 1.2587 + 0.94722Ar -|- 0.2751(r' - i' - 0.5) (A3) 

r'{V, g' - r') = 0.27671 + 0.96792H - 0.22833(g' - r' - 1) (A4) 

We adopted as a measure of uncertainty the standard deviation of the 
residuals of points lying within ±0.5 photographic magnitudes of SDSS 
J0159+0033 (and with a difference between observed and predicted SDSS 
magnitude < 3 mag). Finally, we assume that g — r = 1.0 and r — i = 0.5 
(AB mag) when converting the photographic photometry of J0159+0033 
(which means that the colour terms disappear for SDSS J0159-H0033). Ta- 
ble lAll lists the derived photometric points. 

This paper has been typeset from a T^/ file prepared by the 
author. 
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Table Al. Conversion of DSS photographic plates magnitudes to SDSS photometric system for optical imaging data of SDSS J0159+0033 taken before 1998. 
The source for the observation dates for the photographic plates was taken from https : //archive . stsci . edu/cgi-bin/dss^^late.f inder. 


Epoch 

GSC Band 

Mag 

SDSS band 

Mag (AB) 

1983-09-09 

V 

19.52 ±0.79 

r' 

19.17 ±0.36 

1991-11-02 

N 

18.66 ±0.47 

i' 

18.93 ±0.27 

1992-10-20 

j 

19.94 ±0.36 

9' 

19.70 ±0.22 

1996-10-14 

F 

18.78 ±0.40 

r' 

18.91 ±0.26 
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